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Abstract. We prove semi-classical estimates on moments of eigenvalues of the 
Aharonov-Bohm operator in bounded two-dimensional domains. Moreover, 
we present a counterexample to the generalized diamagnetic inequality which 
was proposed by Erdos, Loss and Vougalter. Numerical studies complement 
these results. 



1. Introduction 

We shall study inequalities for the eigenvalues of the Aharonov-Bohm operator 

(1.1) i/i^ (D - aAo)2 in L^i^). 

Here C is a bounded domain, D :— — iV and aAo{x) :— a\x\^^{—X2,xi)'^ 
is a vector potential generating an Aharonov-Bohm magnetic field with flux a 
through the origin. We shall assume that this point belongs to the interior of the 
simply-connected hull of Q and that a ^ Z, for otherwise aAo can be gauged away. 
On the boundary of fl we impose Dirichlet boundary conditions. More precisely, the 
operator (|1.1[1 is defined through the closure of the quadratic form || (D — q;Ao)u|P 
onCo-(17\{0}). 

Before stating our main results we would like to recall some well-known semi- 
classical spectral asymptotics and estimates for the Dirichlet Laplacian — and 
its magnetic version (D — A)^, A being an arbitrary vector potential. If fl is 
bounded then the spectrum of — A^^ is discrete, and by a classical result due to 
Weyl (see, e.g., jReSj ) one has, as A — + cx3, 

(1.2) tr(-A"-A)I^-i^ // (\^\^-A)ldxd^ = —^\n\A''+' 

i^-^rJJnxR^ 47r(7 + l) 

for all 7 > 0. Note that the right-hand side involves the symbol on the phase 
space 51 X K^. The asymptotics (ll.2p are accompanied by the estimate 

(1.3) tr(-A^^-A)I<i?,-i^ // {\e-Ar-dxd^, 7 > 0, 
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with a universal constant independent of SI and A, and one is interested in 
the sharp value of this constant R-y. In view of (jl.2p the sharp constant obviously 
cannot be smaller than 1, and by an argument of Aizenman and Lieb [AL] , it is 
a non-increasing function of 7. Polya [F] proved the estimate (|1.3p for 7 = with 
constant 1 under the additional assumption that is a tiling domain. His famous 
conjecture that this is true for arbitrary domains is still unproved. Berezin [Blj 
and independently Li and Yau |LiY| (see also [La| ) proved p.3p for 7 > 1 with the 
sharp constant Rj = 1. This also yields the so far best known bound on the sharp 
constant for 7 — 0, namely Rq < 2. Indeed, 

(1.4) tr(-A"-A)0 <(Ai-A)-ltr(-A"-/i)-<(/i-A)-l^|^7|M^ M > A, 

and the claim follows by optimization with respect to fj,. We note that the estimate 
(|1.3p for 7 = and 7 = 1 is closely related to the estimates 

N 

Xn > Po^T^M^^N and "^Xf > pi2n\n\-^N^ 

i=i 

for the eigenvalues of the operator — A^. The form ()1.3p . however, shows the 
close connection with the Lieb-Thirring inequality, see [LTj and also the review 
article |LaW2| . 

We now turn to the 'magnetic' analog of (|1.3p . i.e., where — A^ is replaced by 
the Dirichlet realization of the operator (D — A)^^ in L2(ri) and A is a (sufficiently 
regular) magnetic vector potential. Note that the value of the right-hand side in 
()1.3p does not change if is replaced by ^ — A{x). Hence one is interested in the 
estimate 

(1.5) tr((D - A)2 - A)I < R^'j:^ 1 1 - ^V- dx d^, 7 > 0, 

(27r)^ JJnxR^ 

with a universal constant i?™'*^ independent of fi, A and A. It is a consequence 
of the sharp Lieb-Thirring inequality by Laptev and Weidl [LaWl] that i?™''^ = 1 
for 7 > §• Not much is known about p.5p in the case 7 < |. The Laptev- Weidl 
result and an argument similar to p.4p yield the (probably non-sharp) estimate 
^mag < (|)3/2(-^/(^ _^ ^yp f^^. ^ ^ |^ p^^. ^ = Q and 7 = 1 in particular one finds 

the values 2.1517 and 1.0758, respectively. In [ELoVj the estimate (II. 5p is shown 
to hold for 7 > 1 with constant 1 in the special case of a homogeneous magnetic 
field 

(1.6) A{x) = ^i-x,,x,f. 

It was recently shown in [FLoWj that (|1.5p does not hold with constant 1 if 
< 7 < 1, not even when is tiling. Moreover, the authors determined the 
optimal constant such that (jl.Sp holds for all A under the constraint that A is 
given by (|1.6p and O is tiling. 
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In this paper we shaU consider A corresponding to an Aharonov-Bohm mag- 
netic field and we shall prove the estimate 

(1.7) tTiH^ - A)l < C,{a)-^ [[ {\e-A)ldxd^, 7 > 1, 

with a constant Cj{a) given explicitly in terms of Bessel functions. Even though 
our bound is probably not sharp, it improves upon the previously known estimates. 
Indeed, numerical evaluation of our constant shows that (|1.7p holds for all a with 
constants Co(a) = 1.0540 and Ci{a) — 1.0224 if 7 = and 1, respectively, see 
Section HI We complement our analytical results with a numerical study of the 
eigenvalue of the operator p.ip for five domains: a disc, a square and three different 
annuli. In all cases the estimate (|1.7p seems to be valid with constant 1. We 
refer to Section [5] and Figures [THTOl for a detailed account of the outcome of our 
experiments. 

For the proof of our eigenvalue estimate we proceed similarly as in |ELoVj . 
Indeed, by the Berezin-Lieb inequality, (|1.7p is an immediate consequence of the 
generalized diamagnetic inequality 
(1.8) 

trxn(i?a - A)Ixn < i?7(a) tr xo(-A - K)''_xn for all bounded VL C M^. 

Here Ha '■= denotes the Aharonov-Bohm operator in the whole space. In 
[ELoV] an analogous estimate was proved in the case ()1.6|) with constant 1 when 
7 > 1. The authors conjectured that such an inequality is not true for an arbitrary 
magnetic field, but their counterexample contains a gap; the condition 2) on p. 905 
can not be satisfied by a non-trivial radial vector field, as it was first pointed out by 
M. Solomyak. This gap can be removed by a minor change in the argument, since 
the assumption of radial symmetry is not essential in the proof |ELoj . Nevertheless, 
we feel that our example is of independent interest and sheds some light on the 
particularities of the Aharonov-Bohm operator. What we prove is that the sharp 
constant in (|1.8p is strictly greater than unity, sec Theorem 13.11 We establish this 
by a thorough study of the local spectral density of the operator Ha, see Section [5J 

We mention in closing the papers |BaEvLe] , [EkF] , [H] , [MORj , where Lieb- 
Thirring estimates for the Schrodinger operator Ha + V were obtained. Our esti- 
mates can be seen as a refinement of these estimates in the special case where the 
potential V is equal to a negative constant —A inside and equal to infinity outside 
a bounded domain C M^. 

Acknowledgements. The authors would like to thank A. Laptev for the setting 
of the problem and helpful remarks. The first author is grateful to E. H. Lieb and 
R. Seiringer for their hospitality at Princeton University and thanks them, H. Kalf 
and M. Loss for fruitful discussions. 
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2. The Aharonov-Bohm operator in the whole space 
2.1. Diagonalization 

In this section we recall some well-known facts about the Aharonov-Bohm operator 
in the whole space, see, e.g., |AhBo| . [Ru| . We denote by Ha the self-adjoint 
operator in L2{M.'^) associated with the closure of the quadratic form 

/ |(D-aAo)Mpdx, ite Co°°(M^\{0}). 

Here D := — iV, Ao(x) :— \x\~'^{—X2,xi)'^ and a G K. Moreover, denotes, 
as usual, the Bessel function of the first kind of order i/, see jAbStj . With polar 
coordinates x = \x\{cos9x,sm9x) and similarly for ^, we define 

rtSZ 

and put 

{J^aU)iO-= f Tc,{^,x)u{x)dx, eeK^, 



for u e Co°°(M2). Note that T^{£„x) = (27r)-ie-^«-^ [XbStl 9.1.41], so is the 
ordinary Fourier transform, which diagonalizes Hq = —A. Similarly, one has 

Lemma 2.1. For any a £ R, Ta extends to a unitary operator in L2(R^) and 
diagonalizes Ha, i.e., 

for any u g L2(R^) and f e ioo(R)- 

We sketch a proof of this assertion for the sake of completeness. 
Proof. The orthogonal decomposition 

= 0^„, := {\x\-^/'9{\x\)e^^'- : g G L2(R+)}, 

riGZ 

reduces H^. The part of Ha in i3„ is unitarily equivalent to the operator 

h\n-c,\:=-^ + - ^2 in^2(R+), 

which is defined as the Friedrichs extension of the corresponding differential ex- 
pression on C5"(R+). (We emphasize that in our notation R_|_ means the open 
interval (0, cx)).) The operator 



{<^„g){k) -.^ / VrkJ^{rk)g{r)dr, fc G M+, 

initially defined for g e C^(R-(-), extends to a unitary operator in i2(R+) and 
diagonalizes /i^, i.e., 

{^.f{K)g)ik) = f{k'){<P,g){k), k e M+, 
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for any g G L2{K.+ ) and / G Loo(lli) tTii Ch. VIII]. The assertion of the lemma is 
a simple consequence of these facts. □ 

The proof of the preceding lemma shows in particular that the operators 
Ha and Ha+m with m e Z are unitarily equivalent via multiplication by 6*™^°= 
(a gauge transformation). Hence, without loss of generality we shall assume that 
< a < 1. 

Lemma |2 . 1 1 implies that J {Ha), at least formally, is an integral operator with 
integral kernel 



f{Ha){x,v)= / Ta{Lx)fm'')Ta{U)di 

/(fc2) J|„_„| (fc|x|) J|„_„| (fc|2/|)e"'(^--««)fc dft 



-tT. / /(A)J|„-ai(\/A|:E|)J|„_„|(yA|y|)e™(«--^«)dA. 



On the diagonal this is 

(2.1) f^Ha)ix,x)^— f{X)paiVX\x\)dX, 



where 

(2.2) PoW i>0. 

riGZ 

In particular, Pa{V^\x\) is the local spectral density at energy A. In the following 
subsection we collect some basic information about this function, and in Subsec- 



tions 12.31 and 12.51 we prove some results about the precise asymptotic behavior of 
(12. ip as |x| ^ oo for special choices of functions /. This will allow us to prove that 
the generalized diamagnetic inequality is violated. 

2.2. The spectral density 

Our results are based on a detailed study of the function pa from (|2.2p . We note 
that if a = then po = 1 by [AbSti 9.1.76]. An expression in terms of elementary 
functions is also available |AbSt| 5.2.15] for a = 1/2, 

2 /"^'sins 



(2.3) pi/2(t) = - / ds. 

1^ Jo s 

As t — > oo, Pi/2{t) tends to 1 in an oscillating manner. As we shall see, this behavior 
appears for all < a < 1. The starting point of our study of pa with non-trivial 
flux a is the following 

Lemma 2.2. For any < a < 1, pa is a smooth function on with Pa(0) = 0, 
Pa{t) i as t oo and 

(2.4) p'ait) = Ja{t)Ja-l{t) + Ji_a (t) J_„ (t) , t > 0. 
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Proof. By (LlH 11.2(10)] and [XbStl 11.4.42] one has, for all t > 0, 

Pa{t) = / {Ja{s)Ja-l{s) + Jl-ais)J~a{s)) ds 
JQ 

/oo 
{Ja{s)Ja-l{s) + Ji-a{s)J-a{s)) ds, 

which implies the assertion. □ 

Our next result will not be needed in the sequel, but it helps to clarify the 
behavior of pa and demonstrates the methods which we shall use later on. 

Lemma 2.3. Let < a < 1. As t ^ oo, 

, . sin an cos2t o, 

Pa{t)^l —+0{t-^). 

TT t 

Proof. The asymptotics jAbSt) 9.2.5] 

Mt) . ^- cos (t----) ^ sm (t - 4 - ^) + O(t-) 



the formula (|2.4p and elementary manipulations show that 

, , , 2 / . . „ (2q — 1)^ sinQTT cos2i 

= — sma7rsin2t + '- ^0[r^] 

nt \ At 

Using that Pa{t) = 1 — p'ai^) ds by Lemma [2.21 we obtain the assertion by 
repeated integration by parts. □ 

2.3. Moments of the spectral density 

For any 7 > —1 let us define 

(2.5) cra,j{r) := f (1 - Po,{VXr)dX, r > 0. 



We are interested in the asymptotic behavior of this quantity or, more precisely, 
in the way it approaches its limit. 

Theorem 2.4. Let < a < 1 and 7 > — 1. As r 00, 

1 T./ sinaTT sin(2r - ■i77r) ^-.^^s 



For the proof we note that by Lemma 12.21 and dominated convergence 
lim CTa.^(r) = / (1 - X)^dX ^ 



r — >oo 

In view of 



7 + 1 



cra,j{r)^—-- <j'^^^{s)ds, 
Theorem 12.41 follows via integration by parts from 
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Proposition 2.5. Let < a < 1 and 7 > — 1. As r 00, 
(2.6) 

/ , N r(7 + l)sina7r / / 77r\ di / 77r\ oA 

= ;,+2 , (2r - \) ~ ^ sin (2r - ^) + ^(r-^) j 

w/iere di := i(l - 4((7 + |)2 + (1 - 2a)^). 

We defer the rather technical proof of this proposition to the following sub- 
section. We would like to point out that (|2.6p can be proved in an elementary way 
if a = ^ due to the formula ()2.3p . Indeed, if we simplify further by taking 7 = 1, 
then 

a' , {r) = —[ (1 - A) sin(2\/Ar) dA = — / fc(l - A:^) sin 2fcr dfc 

by integration by parts. This is (|2.6p in this special case. 

2.4. Proof of Proposition 12.51 

We shall need 

Lemma 2.6. Let v > —1/2 and 7 > — 1. Then 

Ju+i/2{kr)J^^il2{kr)k^{l - k^f dk = [ J2^(2fcr)fc(l - k^y+^'^dk 
with 

(0 7\ 1 r(7 + i) 

Proof. By virtue of the series representations |AbSt|, 9.1.10, 9.1.14] for a > 0, 



^ n!r(a + n) V 2 

n— ^ 

/ty-^- (_l)nr(2a + 2n) /t^" 

JairjJa-i W - J „!r(a + n)r(a + n + l)r(2a + n) ^ / 



n=0 

the statement is equivalent to 

(2 8) {''-T'" V (-l)"(^/2)^"r(2. + 2n) ^2.+2„+i^ _ ^2^^^ 

r^"^-! r(7 + l) ^ (-1)'V2" 



y_LiLJ_/ ^2.+2„(^_ ^2)^ + 1/2^^ 

n!r(2j. + n) io ^ ^ 



^ r(7+|),t^on!r(2i. + n)yo 
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an equality which actually holds termwise. This follows by the beta function iden- 
tity 



r(«)r(/3) 



r(a + (3) 

and the duplication formula ^T{2a) = 2'^'^-^T{a)T (a + \) [XbStl 6.1.18]. 
Lemma 2.7. Let /? > and —2 < v <2. The following decomposition holds: 

Mkr)k{l - ef-Hk = /^,,(r) + gp.y{r), 



di . 



□ 



(2.9) 

where g/3,-^{r) = —gp^ir) and 

r(/3) /2\'^+^/^ 



with 
(2.10) 



2^/^ 



cos(r — ro) H sin(r — tq) + 0{r ) 

r 



oo, 



'^0 := TT /3 + + o 



and di 



1 + 



The function has a power-like, non-oscillatory behavior at infinity and 
dominates //j,^ when /3 > |. Its odd parity with respect to u is, however, the 
property vital to us, since it leads to a useful cancellation. 



Proof. The left-hand side in (12. 9|) can be expressed as a generalized hypergeometric 
function. Recall that for p,q E Nq, p < q and ai, . . . , Ofp, . . . , /3g G M \ (— Nq), 



ai,a2 



t'-^At ■ ■ ■ lOtp 
Pi, ■ ■ ■ T Pq 



E 



(Ql)n(«2)n ■ ■ • {Ctp)n 

(/3i)„(/32)n---(/3q)„ n! 



z e C, 



where (a)„ = r(a + n)/r(Q;) = Q;(a + 1) • • • (a + n - 1). By [Lil 13.3.2(10)] 



J^{kr)k{\-k^f-^dk= / J^(r 



7r/2 



' sm f cos 



2/3-1 



r^V{P)V{vl2 + 1) 



-1-F2 



1//2-I-1 
l,/3-l-i//2 + 1 



OsinOde 
1 



2''+ir(i/-hi)r(/3-f i//2 + i) 

Next we use the asymptotics of the generalized hypergeometric function |Lul 
1.3.3(5, 7, 8, 13)], 



T{v/2 + \) 



T{u + l)T{l3 + v/2 + l) 



1F2 



(2/r)'3+''+V2 
v{2/rY+ ^ 

^ 2r(/3) 



J//2 + 1 
+ l,/? + i^/2 + 1 



^d„(*r)-" + e-^('^-'-«)^d„(- 

n=0 



n=0 



-3^0 



1 + 1//2, 1 - 1^/2, 1 - /3 



00, 
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with ro, di as in (|2.10p and do — 1- We emphasize that the ^ sign means equahty in 
the sense of asymptotic expansions, and that 3F0 is not a weh-defined function but 
denotes an asymptotic expansion. We define ly^^r^gjs^^ir) as a finite approximation 
to 3F0. Namely, let K denote the smallest non- negative integer such that 2K > 



/? — I and put 



K 



913 



This function is antisymmetric in v and 



(-1)' 



Mkr)k{l - k^-^dk - gp,,{r) 



2^+1 



m 

'2^ 



2^ 

0+1/2 



ir-ra) ( do + ^) + e^^'-'^") ( do ~ ^ 

\ ir J \ ir 



Or 



-2(K+2) 



do cos(r - ro) + — sin(r - ro) + 0{r ^) 
r 



□ 



as claimed. 

Now everything is in place for the 
Proof of Proposition \2.^ By Lemmas 12.21 and 

<,{r) - 2 
= 2 
= 2c^ 

with from (|2.7p . Now we apply Lemma [HT] with /3 = 7 + | and v ~ ±{2a — 1). 
Using the antisymmetry of gfj^u we find that 

'^a,i{^) — (/7+3/2,2Q-l(r) + /-,.+3/2,l-2a 

and the assertion follows after elementary manipulations from the asymptotics of 
//3,i/ given in Lemma 12.71 □ 

2.5. The Laplace transform of the spectral density 

In this subsection we are interested in the quantity 



k\l ~ k^)'^ p'^{kr)dk 

{Ja{kr)Ja-i{kr) + Ji_„(fcr) J_„(/cr))fc2(l - k^y dk 
{J2a-i{2kr) + Ji_2a(2fcr))A:(l - k^y+^'^dk 



o-Q,oo(r) 



e~^Pa{^/\r)d\ r>0. 



Note that this is essentially the Laplace transform of the function Pa{\f )- 
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Theorem 2.8. For all < a < \ the function la^ao is strictly increasing from to 
1 on [0,oo). In particular, CTa_ao{f) < 1 for all r > 0. 

This theorem shows that the oscillations we observed for aa^-y are no longer 
present. 

Proof. Since pa = Pi-a and hence Ua^oo — fi-a^oo, it suffices to treat the case 
< a < 1/2. By the properties of pa (see Lemma and dominated convergence 
we get (Tq,oo(0) — and <7a,oo(j') ^ 1 as r — oo. Again by Lemma 

/"OO 

Jo 

= 2 fc2g-fc^ (J„(fcr) J„_i(fcr) + Ji_„(fcr) J_„(A:r)) dk, 



and hence according to [Lu, 13.4.1(10)] 
21-2" 



O'a.oo(^) = r 



l-2a 



iFi 



-1^1 



2i-2"r(l-a)' 'V 2 -2a 
In the special case a = 1/2 we note that 

If < a < i we can apply the Kummer transformations [AbSt|, 13.1.27] (note 
that iFi(a, 6; z) = M(a, 6, z) in |AbSt]) to get 



i-2"e-'" 



r2(i-2")r(a) 



li^i 



2a 



1 



'r(i - a 



2- 2a 



By elementary properties of the gamma function this can be rewritten as 

1 



, 2sina7r j_2„ / 1 r, r, 2 

<,oo(0 = e [/( -^a,2-2a,r2 



where [XbStl 13.1.3] 



C/( i-a,2-2a,r2 



1 



3in 2a7r 7.2(1-20) 



1 



r(i-a)r(2a)' V 2a 



2(l-2a) 



-li^l ' 2 



r(a- i)r(2-2a) V 2 - 2a 

U is positive by the integral representation |AbSt| 13.2.5], and this proves the 
theorem. □ 
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3. Counterexample to the generalized diamagnetic inequality 

Following ^ELoVj we consider the question: Which non-negative convex functions 
if vanishing at infinity satisfy 

(3.1) tr xn'piHa)x^i ^ tr xo<i5(— A)xj2 for all bounded domains 57 C M'^ ? 
By (|2.ip the statement (13.11) is equivalent to the pointwise inequality 

/•oo 

ip{X}pa{VXr)dX< (p{X)dX, for aU r e [0,oo). 



(3.2) 

Jo Jo 

Note that (|3.ip is true for the family of functions ip{X) = e~*'^, i > 0; we shall prove 
this (even with strict inequality) in Remark l3 . 2l below. Alternatively, it follows from 
the 'ordinary' diamagnetic inequality (see, e.g., [HuS| ). 



\exp{—tHa)u\ < exp(— 1(— A))|u| a.e., 
and [S Thm. 2.13], since 



u € L2(K^), 



Xnexp ( --^Ha 



= tixncxp{-tHa)xn- 



We point out that the validity of (|3.ip for ip{X) 
that it holds for any function of the form 



-t\ 



, t > 0, actually implies 



(3.3) 



^(A) = 



'{t) dt, 



w> 0. 



In connection with a Berezin-Li-Yau-type inequality one is particularly inter- 
ested in the functions 

(3.4) <^(A)==(A-A)I, 7>1, A>0. 

Note that these functions cannot be expressed in the form (|3.3p . 



Theorem 3.1. Let < a < 1 and let Lp he given by (j3.4p for some 7 > 1, A > 0. 
Then the generalized diamagnetic inequality (|3.ip is violated. More precisely, there 
exist constants Ci,C2 > (depending on a and 7, but not on A) such that for all 
\x\ > CiA-1/2^ 



(3.5) 



ip{Ho,){x, x) - (^(-A)(a;, x) + Aa^^{K) 



sin(2\/A|a;| - ^777) 



|2+7 



A(t-1)/2 
I. T 1 3+7 



With Aa,^{K) := (27r)-2A^/2r(7 + 1) si 



sm an. 



Proof. By (|2.ip and the scaling A 1-^ AA 

A7 + 1 



(3.6) ip{Ha)ix,x)^ 



47r 



(1 - Xr Po.{VAX\x\) dX = ^a,,^(VA|x|), 



where we used the notation ([TS]) . Note that ip{-A){x,x) = (47r(7 + l))-iA'^+^ 
The expansion (13. 5p is thus a consequence of Theorem 12.41 To prove that the 
generalized diamagnetic inequality (j3.1|) is violated, one can consider the domains 
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ri„ := {a; e : \VA\x\ - r„\ < e}, n £ N, with r„ := n{n + ^(7 - 1)) and 
sufficiently small but fixed e > 0. It follows easily from p.Sp that p.ip is violated 
for all large n. □ 

Remark 3.2. The analogous statement (with the same proof) is valid for —1 < 
7 < 1. Moreover, for exponential functions Theorem 12.81 implies that 



1 f°° I 1 ( \x\ 

is strictly less than (47ri)^-'^ = cxp{—t{—A)){x,x). 

The following substitute for p.ip will be useful later on. 

Proposition 3.3. Let < a < 1 and let (p be given by (|3.4p /or some 7 > — 1, 
A > 0. Then for all open sets 17 C 

(3.7) trxnip{Ha)xn < R-t{a) tr xn'p{-^)xn 
with 

(3.8) R^{a) := {j + l)sup f (1 - A)''p„(%/Ar) dA. 



r>0 Jo 

Indeed, this is an immediate consequence of p.6p . This shows as well that 
the right-hand side of (|3.8p yields the sharp constant in p.7p . 

Remark 3.4. The constant i?-y(Q!) is strictly decreasing with respect to 7. Indeed, 
following [ALJ we write, for 7 > 7' > — 1 and < A < 1, 

S(7-7',7' + l)(l-A)^- /' \l-\-tir'fi''-^'-Uf, 



and find, for any r > 0, 

5(7-7',7' + l) / (l-A)>„(VAr)dA 



/" ^(l-A-Ai)'''pa(VAr)rfA^ /i'^-'^'-^d/^ 



where u;(s) := {1~ X)''' pa{V^s) d\. Since ti;(s) < i?^'(a)(7' + l) ^ one concludes 
that 

R^{a) < Ry{a)(j' + + l)B{j - + 1)-^ f\l - a*)^'+V^"^'-' 



Jo 

= i?y (a)(7' + l)-i(7 + 1)B(7 - 7', 7' + l)-'^B{j - 7', 7' + 2) - i?y (a). 
That this inequality is actually strict, follows from the fact that the supremum in 



is attained for some t-q S (0,cxd) (see Theorem 12. 4p and that w{roy/l — fi) is 
non-constant with respect to fi. 



Eigenvalue estimates for the Aharonov-Bohm operator in a domain 



13 



We close this section by giving numerical values for R^{a). Note that 
Rj(a) =2(7 + 1) sup / {\ ~ k^yi pa{kr)k dk 

r>0 Jo 

= supr- / {1 - k^y+^ p'^{kr) dk 

r>0 JQ 

= supr / (1 - k^V+^ {Ja{kr)Ja-i{kr) + .h^„{kr)J^^{kr)) dk. 

r>0 Jo 

The integral can be evaluated by some quadrature algorithm and by Theorem [231 
the supremum is attained for finite r. This allows us to compute approximate 
values of R^{a) (= i?^(l — q)), 





i?-^(0.1) 


i?^(0.2) 


i?^(0.3) 


i?^(0.4) 


i?T,(0.5) 


7 = 





1.01682 


1.03262 


1.04422 


1.05151 


1.05397 


7 = 


1 
2 


1.01027 


1.02050 


1.02781 


1.03241 


1.03395 


7 = 


1 


1.00650 


1.01351 


1.01833 


1.02138 


1.02238 


7 = 


3 
2 


1.00417 


1.00920 


1.01250 


1.01457 


1.01524 


7 = 


2 


1.00267 


1.00642 


1.00874 


1.01019 


1.01065 



In addition to the monotonicity in 7 (see Remark l3.4p the constant R^{a) seems 
to be strictly increasing in a, but we have not been able to prove this. 



4. A magnetic Berezin-Li-Yau inequality 

As in the introduction, let il C be a bounded domain and define in L2(ri) 
through the closure of the quadratic form ||(D — aAo)u|p on C^{fl \ {0}). For 
eigenvalue moments of this operator we shall prove 

Theorem 4.1. Let 0<a<l, 7>1 and SI C be a bounded domain such that 
the operator has discrete spectrum. Then for any A > 0, 

(4.1) tr{H2 - A)I < i?,(a)-^-^|f7|A-+i 



with the constant R^{a) from p.8|) . 
As explained in the introduction, 



^ |f^|A^+i= 1 // {\i\^-K)ld^dx 



47r(7 + 1) ' ' (27r)2 



•■xn 



is the semi-classical approximation for tr(i/^ — A)l. Unfortunately, we can only 
prove (j4.ip with an excess factor R^{a), which is strictly larger than one by The- 
orem 12.41 Based on numerical calculations (see next section) we conjecture that 
(|4?T|) should be valid with R^{a) = 1. 
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Proof. We first note that 

(4.2) tr(iJ^^ - A)_ < tTxnixnHo^xn - A)-Xn- 

Indeed, let {cuj) be an orthonormal basis of eigenfunctions of H^. Then the exten- 
sion ujj of ujj by zero belongs to the form domain of x^H^xn and 

tr(i?^ - A)_ - ^((i/i^ - A)^„^,)_ = Y.^{xnH^Xn - A)i;,>,)- 

j 3 

< ^{{xnHaXn - A)_a;j,a;j) = trxn(xai?aXa - A)_xn- 

3 

Now let if he a convex function of the form ip{X) — /(A — ii)-w{fi) dfj, for some 
w>0. Then it follows from that 

tTifiH^) < trxo¥'(xoi?aXo)xo, 
and hence by the Berezin-Lieb inequality (see |B2j . [L] and also |LaSa| . |Laj ). 

tr<y5(-ffa) < trXJi¥'(xo-ffQXn)xo < i^: Xn'p{Ha)xQ■ 
As already noted in Remark l3.41 the function ip{X) = (A — A)l is of the considered 
form. Noting that in this case tr xn'fii— ^)xn = (47i'(7+l))^"'^|r2|A'''+-'^, the assertion 
follows from Proposition 13.31 □ 



5. Numerical experiments 

5.1. Numerical evaluation of the magnetic Berezin-Li-Yau constant 

In this section we present a numerical study, somewhat in the spirit of [LTi Ap- 
pendix A], of the constant Rj{a) appearing in (|4.ip for on various domains. 
We consider the unit disc, a square and three annuli: 

{x -.Ix] <1}, C := {x : 1< |x| < 1.1}, 

B := {x : max(|xi|, \x2\) < 1}, D := {x : 1 < \x\ < 2}, 

£; := {x : 1< |a:| < 11}. 

To expose possible diamagnetic effects we perform identical experiments for a = 
and a = 0.2 throughout. If the eigenvalues Ai, A2, . . . are known, the task is to 
determine 

(5.1) i?^ = supr^(A), where r^(A) ^^^^ ^ (A - A,r . 

Since the sum is simply the counting function if 7 = 0, tq is decreasing for all 
A but the eigenvalues. In the physically important case 7=1, the quotient ri is 
continuous but r[ has jump increases at the eigenvalues. Figures [T] and [H show, 
respectively, the schematic behaviour of near the bottom of the spectrum and for 
large A in the case where — I (in fact, the spectrum of Hq has been considered). 
Figure [3] depicts ro(A) for an operator with > 1 (see [HI Thm. 3.2]), namely 
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the Schrodinger operator Ha — \x\ ^. The x symbols on the x axis indicate the 
eigenvalue loci. 

5.2. Computation of the eigenvalues 

By separation of variables the spectrum of is the union of the spectra of the 
family of one-dimensional problems 

-u" + ~ ~ ^^^ u ^ \u, w(l) = 0, ueL^{0,l), 

parametrized by n g Z. \/rJ\n-a\i\^''') is an eigenfunction provided A is chosen 
to that J\n-a\ (VA) = 0, i.e., the eigenvalues are squares of the zeros of the Bessel 
function of the first kind. Using standard numerical procedures for this task we 
compute all eigenvalues below 10,000 (approximately 2,500) of Hq and all eigen- 
values below 50,000 (approx. 12,000) of H^2- The lowest eigenvalues of the two 
operators can be viewed in Figure 31 The reader interested in the dependence of 
the eigenvalues on the parameter a may find Figure [H] useful, where the lowest 
eigenvalues of as functions of a G [0, ^] are shown. We recall that and 
Hi_a share the same spectrum. 

Similarly, the eigenvalue problems for H^, and are reduced to 

,, (n — a)^ — 1/4 , , , , , 9, N 

-"^ +- ^ —u = \u, u{ri) = = u{r2), u£L^{ri,r2). 

In order that the general solution •\/r(ci>/|ri-a| (V^?') + C2i^n-a| (n/At")) satisfy the 
Dirichlet conditions we need to fix the ratio ci/c2 and take A such that 

^|n-a|(VAri)l|„_„|(\/Ar2) = J|„-a|(VAr2)Y|„_Q|(VAri). 

Equivalently, A = x^/r^ is an eigenvalue whenever a; is a zero of the cross-product 
J^(t)Y^(lit) — Ji,{pt)Y^{t) for fi = Tijrx^ v — \n — a\. Since the amplitude of the 
oscillation grows approximately as jv^ the computation of the zeros is nontrivial 
for large orders. By searching in the negative t direction and making use of the 
simple fact that the asymptotic spacing of the zeros is 7r/(/i — 1) |AbSt[ 9.5.28], 
we circumvent much of these difficulties in comparison to routines such as Math- 
ematica's BesselJYJYZeros. This allows us to compute, for a = and 0.2, all 
eigenvalues below 10,000 (approx. 400) of , all eigenvalues below 1,000 (approx. 
550) of iJf and all eigenvalues below 300 (approx. 9,000) of . 

Hq is simply the Dirichlet Laplacian on [—1, 1]^, so for all k.l G N, 

, x + l\ . / y + l 

KIT Sm ITT 



is an eigenfunction and \{k'^ + 1'^) an eigenvalue. For a ^ the absence of radial 
symmetry prevents us from solving the eigenvalue problem for exactly. We 
therefore resort to the finite element method as implemented in Comsol Mul- 
TIPHYSICS to compute approximate eigenvalues. A 8,800-element mesh, gradually 
refined near the singularity, is used to compute the 250 lowest eigenvalues, all less 
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than 840. Providing an a priori bound on the computational error would be beyond 
the scope of this article; it is an encouraging fact, however, that the approximate 
eigenvalues of do not differ from the exact ones by more than 0.8 %. Figure[5] 
is the analogue of Figure [3] for . 

5.3. Outcome of the experiments 

We are obviously constrained to computing (jS.ip with A restricted to a finite 
interval but provided this interval is sufficiently large it should be possible to tell 
the behaviors depicted in Figures [2] and [3] apart. Particularly if the function 
tends steadily to 1 from below we can conclude with certainty that its supremum 
will eventually be this number in virtue of the eigenvalue asymptotics. It turns 
out that all the operators we consider, whether a = or 0.2, seem to obey a 
semi- classical eigenvalue estimate with unit constant for both 7 = and 1. 

What however differs is the rate of convergence to 1, e.g., r^ increases slower 
for than for since |C| < \E\, implying that the point spectrum of the 
former operator is sparser. To see the influence of a on we invite the reader to 
compare Figures [S] and [51 Note that r^ has the most irregular behavior in the 
cases a = and a = 0.5 due to the high degree of degeneracy of the eigenvalues, 
cf. FigureEl Moreover we have the impression, from looking at Figures fTI, 171 and fTUl 
that the shape of the domain has a limited influence on r^ . All three plots suggest 
that i?o = 1 but this has not been proved so far even for Hq, the Laplacian on 
the unit disc! 
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Figure 2. r-y(A) for for large A 
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Figure 5. Lowest eigenvalues of Hq (•) and of 2 (x) 
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Figure 6. Lowest eigenvalues of as function of a 





Figure 8. for near the bottom of the spectrum 



